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Abstract 

We review the status of the semiclassical trace formula with emphasis on the particular 
types of singularities that occur in the Gutzwiller-Voros zeta function for bound chaotic 
systems. To understand the problem better we extend the discussion to include various clas- 
sical zeta functions and we contrast properties of axiom-A scattering systems with those 
of typical bound systems. Singularities in classical zeta functions contain topological and 
dynamical information, concerning e.g. anomalous diffusion, phase transitions among gen- 
eralized Lyapunov exponents, power law decay of correlations. Singularities in semiclassical 
zeta functions are artifacts and enters because one neglects some quantum effects when de- 
riving them, typically by making saddle point approximation when the saddle points are 
not enough separated. The discussion is exemplified by the Sinai billiard where intermit- 
tent orbits associated with neutral orbits induce a branch point in the zeta functions. This 
singularity is responsible for a diverging diffusion constant in Lorentz gases with unbounded 
horizon. In the semiclassical case there is interference between neutral orbits and intermit- 
tent orbits. The Gutzwiller-Voros zeta function exhibit a branch point because it does not 
take this effect into account. Another consequence is that individual states, high up in the 
spectrum, cannot be resolved by Berry-Keating technique. 



I. Introduction 



Research in quantum chaos has to a large extent been focused on conceptual questions 
such as how is classical chaos revealed in the h ^ limit of quantum mechanics and what 
are the properties of quantum spectra of classically chaotic systems. During the last decade 
another hot topic has been periodic orbit quantization. Model calculation of quantum spectra 
using periodic orbit techniques are greatly simplified if one restrict oneself to systems with 
two degrees of freedom. Applications to three dimensional systems are then restricted to cases 
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where some internal symmetry reduces the system by one degree of freedom. One example is the 
diamagnetic Kepler problem, see e.g. Q and references therein. The simplest conceivable two 
dimensional systems are billiards. Billiards are of great use in theoretical work on semiclassical 
mechanics. But billiard like systems, of real interest, and of any imaginable shape can now be 
constructed in the laboratory, there are several examples in the present volume. So the stage 
is set for a fruitful interplay between quantum chaos and mesoscopic physics. 

To make a successful semiclassical computation one needs to control the semiclassical limit 
of the Green function. One can then compute anything from conductance in mesoscopic devices 
[|[ to wave functions. Another information contained in the Green function is the quantum 
spectrum, which can be uncovered by taking its trace. In bound systems the Green functin 
can be expressed in terms of the eigenfunctions and eigenvalues of the Hamiltonian 

i * 

and the trace reads 



trG{x,x',E) = J G{x,x,E)dx = ''^— — 

i 

whose semiclassical limit for chaotic system was derived by Gutzwiller 



1 CO -. 

trG{E) = go{E) + -J2TpJ2 Ms,/n^>.,^] (3) 

The index p labels the primitive periodic orbits, Sp is the action integral along the orbit, Tp 
(= dSp/dE) its period, Mp is the linearized Poincare map around the orbit and fip is a phase 
index. Finally, go{E) provides the mean level distribution, which may be interpreted as the 
contribution from classical orbits of zero length. 

The trace formula can be recasted to the Gutzwiller- Voros zeta function^ 

00 / „45p/?i-/.pf ] \ 

Zgv{E)=\{\{ l-^^nT^ • (4) 

The oscillating part of the trace formula is given by the logarithmic derivative of the zeta 
function. The only continous time systems we are going to consider are billiards. The action 
Sp is then related to the geometric length Lp of the orbit as Sp = Lpk where k is the momentum. 
Usually we are going to consider zeta functions as functions of the complex variable s = ik. 
This paper is devoted to the study of classical and semiclassical zeta functions. In particular, 
we will be interested in the presence of various kinds of singularities. This might seem a highly 
technical problem but by putting this theme in focus we will shine some new light on many 
of the problems that are associated with semiclassical quantization. The paper is divided into 
two parts. First a nonechnical review of the subject. Then follows a series of more technical 
discussions to illustrate various aspects of the preceding story. 
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II. A Story - The rise and fall 
of periodic orbit theories in quantum mechanics 



Our story will begin with Berry's 'tour de force' paper Quantizing a Classically Ergodic 
System, Sinai's Billiard and the KKR method. Centered around a numerical method Berry 
discusses virtually all issues now known under the phrase of Quantum Chaos, like e.g. level 
statistics and periodic orbit expansions. Taking the semiclassical limit of level density he 
obtains the Gutzwiller trace formula for the unstable orbits and the Berry- Tabor like expression 
for the neutral (or marginally stable) orbits. 

At the time it was premature to start from the trace formula and compute (semiclassical 
approximations of) the eigenvalues. Berry's numerical scheme was successful because of a major 
trick; a transformation of the badly converging sum over lattice points to a highly convergent 
sum over points in the dual lattice. One of the major problem with periodic orbit expressions 
is to make the corresponding trick for the Gutzwiller sum. We will discuss this problem from 
several points of view in the following. 

Later, Gaspard and Rice were inspired by Berry's method to construct a numerical 
scheme for determining scattering resonances in the 3-disk scatterer, and similar systems. 
Prom one point of view this is an easier problem (from another, more conceptual it is not), the 
disk lattice is finite and there is no need to go over to the dual lattice. 

The 3-disk is easy to deal with also from a semiclassical point of view, as was realized 
by Cvitanovic and Eckhardt Q. The Gutzwiller- Voros zeta function^ closely resembles zeta 
functions studied in classical problems. The 3-disk problem is ideal for application of classical 
zeta functions. It has a simple topology, provided the disks are enough separated, the dynamics 
may be mapped to a dynamics among a finite set of symbols. This is known as symbolic 
dynamics. The open 3-disk is also hyperbolic; the stabilities of periodic orbits are exponentially 
bounded with their length. A system with these two properties is said to be of Axiom-A type. 
This property gives nice convergence to the Dirichlet series obtained by expanding the Euler 
product representation of the zeta function (^) , this is usually called a cycle expansion ||8|, 0] . 
The classical resonances as well as the quantum resonances can thus be efficiently computed 
by a moderate number of periodic orbits. 

Still, the situation is not ideal. The Gutzwiller- Voros zeta function does exhibit poles |^] 
which limits convergence of cycle expansions far down in the complex fc-plane. A theorem 
of Rugh[10| says that a certain type of zeta functions, Fredholm determinants of evolution 
operators with multiplicative weights, are entire for Axiom-A systems. A Fredholm determinant 
may look something like 

nz) = n n (1 - ^^r^^ (5) 

P m=0 

The Gutzwiller- Voros zeta function lacks the right structure among the higher m factors to 
make it entire. There is a lot of interesting work [11, 12| on improving the Gutzwiller- Voros 
zeta function by manipulating these higher factors. The by far most interesting is the quasi 
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classical theory of Cvitanovic and Vattay |12]. Starting from the hydrodynamical formula- 
tion of Quantum Mechanics of Madelung they construct a multiphcative evolution operator 
whose Predholm determinant should contain the semiclassical eigenvalues in its spectrum. The 
derivation avoids the multitude of stationary phase approximations in the standard derivation 
of the trace formula. Still, the quasiclassical Fredholm determinant contains extra zeros which 
have nothing to do with the quantum spectrum, and they seriously slow down convergence. So 
the surprising thing is that, if interpreted as an asymptotic series, the Gutzwiller-Voros zeta 
function is in many respects still superior to its competitors |13|. 

Due to the simplicity of the 3-disk scatterer it is an ideal testing ground for various 
corrections to the Gutzwiller-Voros zeta function, like diffraction-corrections [14, 15| and h 
corrections 1 16]. Note that the circumstances that make the 3-disk classically and semiclassi- 
cally nice also make diffraction corrections small; semiclassical formulas fail if the point particle 
scatter in an extreme forward angle, or if it passes extremely close to a disk and this is never 
the case in the (enough separated) 3-disk. The main diffraction effect for this system is then 
due to creeping and this is indeed a small effect. 

When including diffraction effects in axiom-A systems one can say that the classical phase 
space is enlarged to include diffractive orbits associated with some discontinuity such as a point 
or wedge or circular disk. This can be done without risking that stationary points come to 
close. We will see that this is very different from the problem intermittency is going to present 
us with in bound systems. 

We have dwelt on Axiom-A scattering systems at some length, mostly to appreciate how 
different from bound systems they are. When we now are turning our interest towards bound 
systems we will encounter singularities that have nothing to do with the higher ni factors in 
Fredholm determinants, so we will frequently omit them. Their subtle role in bound systems 
is far beyond the scope of this paper. 

If the open 3-disk system was an archetype for open chaotic scattering systems the Sinai 
billiard is a suitable archetype for bound chaotic systems. It is in many ways generic (see below) 
but is simple enough to allow some analysis. This system is either represented as a circular 
scatterer inside a square or as an infinite lattice of disks, cf. fig 1. The latter representation 
is an example of a Lorentz gas. Lorentz gases are among the simplest conceivable systems 
exhibiting statistical behaviour such as diffusion p^ . 

Bounded system are generally not of Axiom-A type. We will now discuss various deviations 
and how they affect the analyticity of the zeta functions. 

The topology is generally extremely complicated. Given a bound system, there is a priori 
no reason to believe that there should exist a simple symbolic dynamics. Often one can find a 
symbolic coding but the grammar rules are complex, in the sense that the number of forbidden 
substrings increases exponentially with length. Mathematically, such a symbolic dynamics is 
known as an infinite subshift. The basic mechanisms for pruning of the symbolic dynamics can 
be very different. In dispersive billiards pruning enters because a trajectory from one disk to 
another can be obstructed by an intermediate disk. 
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Very little is known about how pruning affects the analyticity of zeta functions. In ref, 
we propose that the lack of finite symbolic dynamics may induce a natural boundary in 
the complex plane through which analytic continuation is not (even in principle) possible. We 
discuss this in little more detail in section Ilia. 

Intermittency is another abundant property of bound chaotic systems. We will consider a 
strong form of intermittency associated with the presence of neutral orbits in the system. Orbits 
accumulating towards this neutral orbit will have Lyapunov exponents tending to zero. Many 
ergodic billiards, such as the stadium|19] and the Sinai billiard have this property. The effect 
is relevant for Hamiltonian systems with mixed phase space because the boundaries of stable 
islands are neutrally stable and will be surrounded by intermittent orbits. So the neutral orbit 
in a chaotic billiard may be looked upon as a stable island of zero area [|l9|| . Strong intermittency 



appears to imply a branch point singularity 21 1 at the origin (s=0) as discussed in section 
Illb. It is possible to construct bounded billiards without neutral orbits. One example of is 
(the unit cell) of a Lorentz gas on a triangular lattice, provided the disks are large enough, see 
[ p^ . Such systems often exhibit some weaker form of intermittency. 

In classical zeta functions, these singularities play a dual role. They hamper convergence of 
cycle expansions seriously. But they are also important carriers of dynamical and topological 
information such as power law decay of correlations, anomalous diffusion and phase transitions 
among the generalized Lyapunov exponents [p^j, and should not be considered as undesired 
mathematical obstacles. In section IIIc we, as an example, discuss the relation between the 
singularity of the zeta function and anomalous diffusion in the Sinai billiard. 

Singularities will be present in virtually any type of zeta function, and in particular the 
Gutzwiller-Voros zeta function (see section Illb). This zeta function is supposed to approxi- 
mate the spectral determinant A{k) (apart from the smooth Weyl factor). This function obeys 
a functional equation A(k) = A{—k) and the singularities we have discussed violate this func- 
tional equation. In the Gutzwiller-Voros zeta function they are artifacts and enter because 
one neglects some quantum effects when deriving it, like making saddle point approximations 
when the saddle points are not enough separated. This problem is intimately associated with 
intermittency. In section Ille we demonstrate that the problem with interfering stationary 
points actually occur in the Sinai billiard and this is what lets a branch point enter the zeta 
function. 

If the system under study is hyperbolic the situation is somewhat better, but one has 
still to bother about pruning. The main pruning mechanism in dispersive billiards is, as we 
said, due to trajectories between two disks being obstructed by a third disk. Orbits that are 
nearly pruned are therefore greatly affected by diffraction effects so pruning and diffraction goes 
hand in hand and must be dealt with in a coherent way, let's call such a hypothetical scheme 
Quantum pruning. Pruning in time reversible system can often be described by a pruning front 
[ p^ . Intuitively we can say that the quantum counterpart is not sharp as the classical one. 
Recall that we expect that a generic pruning front may induce a natural boundary in classical 
zeta functions, so the fuzziness in the quantum counterpart should prevent such a singularity 
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to occur. 

How serious the problems of the Gutzwiller-Voros zeta functions are depend on the method 
by which one computes the zeros. There are mainly two approaches. The first is the cycle 
expansion approach, previously discussed, used mainly in open Axiom A systems. This relies 
heavily on strong convergence properties of the zeta functions. An alternative for bound systems 
is to use the method of Berry-Keating [ 26 1. The main idea is the following. The existence of a 
functional equation for the spectral determinant imply that one only needs to know the length 
spectrum (the fourier transform of the level density) up to half the Heisenberg length]^ in 
order to compute the spectrum. The most accurate computation of spectra of the Gutzwiller- 
Voros zeta functions for bound systems has been performed using this approach[27, 28|. The 
advantage of the method is that one only needs a finite number of periodic orbits to compute 
a finite number of levels, and one only needs to scan the real energy axis. The approach of 
Berry-Keating assumes that the spectral determinant may be expressed as a Dirichlet series. 

There are now growing evidence p|, ||] that even the Berry-Keating method will eventually 
break down for billiards. In ref. [^] the authors argue that the majority of orbits in the Sinai 
billiard with length L > 0(A;^/'^) is affected by diffraction and that these corrections are of the 
same order as the original contribution. Recalling that the Heisenberg length in billiards scales 
as Lh = 0{k) this implies that the Berry-Keating scheme will break down as /c — > oo, provided 
that diffraction effects cannot be accounted for. The occurrence of a branch point singularity is 
actually just another side of this coin. It indicates that the tail of the Dirichlet series suffers so 
badly from intermittency and the nearby stationary points so as to make the series divergent (a 
branch point at A; = makes the real /c-axis where we want to find our eigenvalues the border 
of convergence). The entanglement between neutral and intermittent orbits implies that the 
representation of the spectral determinant in terms of a Dirichlet series is lost. There is thus 
no obvious implementation of the Berry-Keating method. 

The problem with interfering saddle points has nothing to do with the singular nature 
of billiards. In smooth chaotic potentials it is even harder to realize a complete symbolic 
dynamics. The abundance of pruning is deeply associated with the existence of small stability 
islands. These islands will immediately cause problems with intermittency and close saddle 
points. 

It should be stressed that such a failure of predicting individual eigenstates does not render 
the trace formula useless. One may take two approaches. One is to take the trace formula 
more or less as it stands and ask for the limitations of it. What information can be gained 
from this simple and beautiful formula? To use it as a theoretical tool one also needs a better 
understanding of the asymptotic properties of periodic orbits |37, ^] far beyond the Ozorio- 
Hannay sum rule |3^. The other approach is to try to improve the trace formula in various 
ways. The two approaches are well motivated and are of course not independent of each other. 
It is an intellectual challenge to improve periodic orbit theories to predict eigenstates higher 

^the length corresponding to one mean level spacing 



6 



and higher but it is perhaps not always what apphcations ask for. 

Among all systems exhibiting hard chaos, systems with neutral orbits are among the worst, 
without neutral orbits the situation would probably be better. The conceptual problems of 
systems with neutral orbits are similar to those of generic systems with mixed phase space, 
and are thus of fundamental importance. The relative simplicity of the Sinai billiard make it 
an ideal model system for addressing these issues. 



III. Five Illustrations 

We will now illustrate the story in the previous section with a series of more technical 
discussions. We are trying to understand the sense of singularities in zeta functions in the 
classical case and the nonsense of them in the semiclassical case. We will start our expose by 
considering classical zeta functions for discrete time systems (i.e. maps), move on to continous 
time systems (in particular the Sinai billiard) and finally discuss semiclassical zeta functions. 

a. Symbolic dynamics in bound chaotic systems 

We begin by considering one dimensional maps x i— > /(x). Much dynamic information, like 
different kind of entropies and generalized Lyapunov exponents, is encoded in the following 
one-parameter family of zeta functions |31] 



ymz) = Y[{i--^) . (6) 

The product in (^) runs over all primitive periodic orbits p, having period rip and stability 
Ap = with Xp being any point along p. If /3 = 1 no metric information enters the zeta 

which is now called a topological zeta function 

l/Ctop{z) = ll{l-z^^) . (7) 

p 

The leading zero zq of l/Ctop(-z) (the one with smallest modulus) and the topological entropy 
h are related by /i = — log zq. 

We will restrict our attention to the topological zeta function for unimodal (one-humped) 
maps with one external control parameter f\{x) = \g{x). Symbolic dynamics is introduced by 
mapping a time series . . . Xi-iXiXiJ^i . . . onto a sequence of symbols . . . Sj-iSjSj+i . . . where 

Si = L Xi < Xc (8) 
Si — C Xi — Xq 5 (9) 
Si = R Xi > Xc (10) 

and Xc is the critical point of the map (i.e. maximum of g). The kneading sequence is the 
itenary of the critical point. The allowed symbol sequences can be determined from it. All 
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unimodal maps (obeying some further constraints) with the same kneading sequence have the 
same set of periodic orbits 

If /3 / 1 in (^) the individual details of f{x) is reflected in the zeta function, but for tent 
map 

[ A • (1 — x) X G (1/2, IJ 

(where the parameter A G [0, 2]) the general zeta function (|6|) is obtained from 1/Ctop by simple 
rescahng, 

l/C(/3,^) = l/Ctop(^/Ai-^) . (12) 

The topological entropy is h = — log A. 

The set of periodic points of the tent map is countable. A consequence of this fact is that 
the set of parameter values for which the kneading sequence is periodic or eventually periodic 
(preperiodic) are countable and thus of measure zero and consequently the kneading sequence 
is aperiodic for almost all A. For general unimodal maps the corresponding statement is that 
the kneading sequence is aperiodic for almost all topological entropies. 

For a given periodic kneading sequence of period n, I_x = PC = siS2 ■ ■ ■ Sn-iC there is a 
simple expansion for the topological zeta function. Now let Oj = 1 if Sj = L, and Oj = — 1 if 
Si = R. Then the expanded zeta function is a polynomial of degree n 

n-l 

yCtopiz) = - = (1 - ^) • E , (13) 

P 1=0 

where 

n 

bn = l[a^ . (14) 

i=l 

Aperiodic and prep er iodic kneading sequences is accounted for by simply replacing n hj oo. 
An important consequence of (13) and (|14|) is that the sequence {hi} has a periodic tail if and 



only if the kneading sequence has one (their period may differ by a factor of two though). 

The analytic structure of the function represented by the series ^ biZi depends on whether 
the tail of {bi} is periodic or not. If the period of the tail is we can write 

1/Ctopiz) = p{z) + q{z){l + + . . .) = p{z) + ^44v ' (15) 

for some polynomials p{z) and q{z). The result is a set of poles spread out along the unit circle. 
This applies to the preperiodic case. An aperiodic sequence of coefficients would formally 
correspond to an infinite N and it is natural to assume that the singularities will fill the unit 
circle. There is indeed a theorem ensuring that this in a sense is the case (provided the 6j's can 
only take an finite number of values). The unit circle becomes a natural boundary. A function 
with a natural boundary lacks an analytic continuation outside it. 
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What happens with z) for /3 7^ 1 if we make our unimodal map nonhnear but still 

hyperbolic (that is |/'| > C > 1)? All that is known is that the the zeta function is analytic 
inside a certain radius and an obvious guess is that this radius is limited by a natural boundary 
in the generic case. It is natural to expect that the natural boundary moves outward from the 
unit circle as (3 is decreased, and probably it gets deformed. 



Even less is known if / is intermittent, some results have recently been made |34|. Except 
for the branch cut (see section Illb) the zeta functions (/3 = 0) is holomorphic in a certain 
region enclosing the unit circle. 

If the system is time reversible one expects pruning to be described by a pruning front 
which is the two-dimensional generalization of the kneading sequences [p^]. Moreover, if the 
system is time continous the expanded zeta function is a Dirichlet series rather than a power 
series. Knowing that natural boundaries are abundant also in Dirichlet series we have at least 
to focus the possibility of the existence of natural boundaries in zeta functions of systems with 
no simple symbolic description. 

b. Intermittency 

Intermittency means that the system switches between chaotic behaviour and quasi-integrable 
behaviour. We will consider a strong form of intermittency associated with the presence of (one 
or more) neutrally stable orbits. To illustrate how this affects the analyticity of the zeta func- 
tion it is instructive to consider a one dimensional map x f{x) on the unit interval, with 

^, , [ x + 2'x^+' < X < 1/2 

^(") = \ 2.-1 l/2<x<l ' ^''^ 

where s > 0, see fig 2. For s = the map is just the binary shift map, which is uniformly 
hyperbolic, but for s > it is intermittent; the fix point x = is neutrally stable: /'(O) = 1. 
The map admits a binary coding, as before we associate the letter L with the left leg, and R 
with the right leg. The neutral fix point now corresponds to the periodic orbit L. The symbolic 
dynamics is by construction complete. 

We will consider the zeta function (^ with the neutral fix point pruned 

i/C/(^)=n(i-rvi^) ■ (17) 

We now suppress the parameter /3. The new index indicates what map the zeta function refer 
to. 

A power series representation of l/(" is obtained by expanding this product: 1/Cf{z) = 
J2n(^i^^- The nature of the leading singularity will be reflected in the asymptotics of the 
sequence To get an idea what this asymptotic behaviour may be we consider the zeta 

function of closely related map /, a piecewise linear approximation of /. We define / as a 
continous function, coinciding with / on a sequence of points /(c„) = f{cn) where the c^'s are 
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the inverse images of the critical point (cf fig 2) 



CO = 1/2 (18) 

Cn+l = fE\cn) , (19) 

and hnear in the intervals [c„+i,c„]. 

It is relatively easy (see e.g. [^]) to show that the sequence {c„} has the asymptotic 
behaviour 

Cn ~ n~^/' n^oo . (20) 
The construction of / gives it a very simple cycle expansion 

oo n-\-l 

1/C;(.) = 1-E|= = E«"^" • (21) 
The stabilities Ap^^ are simply related to the c^'s 

A^ = 2 ""~""+^ , (22) 

^ « Co - Ci 

with the asymptotic behaviour 

Al^~^^^+'^/' • (23) 

The asymptotic behaviour of the coefficients are thus Oi = 0(1^(^^+1)/^) and it seems likely that 
the same holds for the sequence {ai}. This suggest that \/C,f{z) contains a singularity of the 
type 

(1 - zT-^ « iv 

(1 - z)"^Mog(l - z) aeiV+ , ' 



with 

(l-/?)(^ + l) 

a{f3,s) = , (25) 

s 

as can be realized through the Tauberian theorems for power series. 

We will now move on to the Sinai billiard. In billiards with continous time one cannot 
study zeta functions in complex z since the lengths Lp are not integer multiples of some unit 
length. Instead we replace z by exp(— s) and formulate the the zeta functions 

^(«) = 11(1 - (26) 

where p runs over all primitive unstable cycles. A semiclassical zeta function is formally 
obtained as /3 = 1/2, we only need to reinsert the Maslov indices (which can be done properly 
by a multiplicative weight). 

The presence of a simple symbolic dynamics in the previous example made it easy for us 
since (at least in the piecewise linear approximation) we only needed to consider one sequence 
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of periodic orbits accumulating at the neutral fixpoint. In the Sinai billiard there are neutral 
orbits in a (finite) number of directions. A sequence of periodic orbits accumulating towards 
a neutral orbit is indicated in fig 1, let us denote the sequence g„. The stability of the cycles 
obey a power law Ag^ ~ -n? asymptotically and since Tq^ ~ n the local Lyapunov exponent 
log Ag^/Tg^ goes to zero and we can again expect a branch cut. But the order of this singularity 
cannot be estimated as easily as in the previous one dimensional map. The reason is the lack of 
simple symbolic dynamics. There is a simple symbolic coding |37] but it is heavily pruned for 



cycles having segments close to the neutral orbits. We must rely on more elaborate methods 
to investigate the singularity. 



In a series of paper pq, p7^, pi], pq, El, 24 1 we have investigated this problem within the 



frame work of the BER (acronym for Baladi-Eckmann-Ruelle) approximation |39|]. The basic 
idea is that for intermittent systems one can usually find a surface of section whose associated 
Poincare map is not intermittent, ideally it is exponentially mixing. An approximate zeta 
function for the flow through this map can be formulated as a an average over the surface of 
section 

^BSK(s) = l-(e-^^^("^)|A(x,)|^),,e,o, , (27) 

where As(xs) is the traveling length to the next intersection with the surface of section and 
A(xs) its associated instability, both being functions of the surface of section coordinate Xg. 

In Hamiltonian chaotic system this average is easily performed since the invariant density 
is uniform. It is convenient to define 

Pp{/^) = {5{/^-/^s{Xs))\K^s)\^)sos . (28) 

For /3 = this is just the (probability) distribution of recurrence times to the surface of section. 
The BER zeta function is obtained by a simple Laplace transform 

POO 

Zber{s) = 1 - / e-^%(A)c?A (29) 
Jo 

In the Sinai billiard the obvious choice of surface of section is the disk. The function pp{^) 
of the disk-to-disk map can be controlled well: for /3 = we get in the limit i? ^ 

Po(A) ~ 

1 ||(2e + e(4-30iog(e) + 4(e-i)2iog(e-i)-(2-o'iog|2-ei) e>i ' 

where ^ = /S.2R. 

For finite R we only have simple expression for the tail 

P/3(A) - ^3-3/?/2 • (31) 
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The prefactor can be computed exactly for any disk radius R but this is not our concern here. 
It is also well known the the mean is given by 



Ap^=o{A)dA = ^ + 0{R) , (32) 



an expression we will need in section IIIc. 

A power law tail of p ~ A~" thus gives a leading singularity 



log s a G A^"*" 



„a—l i„„ „ „. ^ i\T+ 

(33) 



to the BER zeta function. 

How accurate the BER approximation is depend on the underlying map. For the Sinai 
billiard the disk to disk map is hyperbolic. It is not known if it i exponentially mixing but 
the correlation function can at least be bounded by a stretched exponential. We expect 
that the first few terms in the generalized series expansion is exactly reproduced by the BER 



approximation, at least we expect that |24] 



Z{s) = Zber{s)+0{s^) , (34) 



but higher terms are also described well by the BER approximation |37|. For /? = this would 
imply that the leading singularity log s is exactly reproduced by the BER approximation. 
For the semiclassical case /? = 1/2 we get a leading singularity of the form s^/^ and thus an 
algebraic branch point at s = 0. The introduction of Maslov indices will affect the size but not 
the order of the singularity. 



c. The role of singularities in classical zeta functions 
What are classical zeta functions good for? 

The most well known application is the relation between the zeros of l/C/3=o and the decay 
of correlations. This problem is rather complicated mathematically and we will instead consider 
the somewhat simpler problem of computing averages of chaotic systems In particular 

we will be interested in computing the diffusion constants but we'll start by considering more 
general averages. 

We assign a weight w{xQ,t) to the trajectory starting at phase space point xq and evolving 
during time t to point x{xo,t). There will be an important technical restriction on w{xo,t): it 
must be multiplicative along the flow, which means that w{xo, ti+t2) = w{xo,ti)w{x{xo,ti),t2). 
The phase space average of w{xo,t) may be expanded in terms of periodic orbits as |^3[ 

limKxo,t))= hm^r,^^; \ , (35) 

p r=l I f' 
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where r is the number of repetitions of primitive orbit p and Wp is the weight integrated along 
cycle p. This trace may be formulated in terms of zeta functions in the following way 



hm {w{xo,t)) = lim — / e^'^ffids , (36) 



with the zeta function 

^«.(s)=n(l-^P^j • (37) 

In the previous section we studied the thermodynamic weight Wp = |Ap|^ but we will now 
turn to the problem of diffusion of interest for the present volume. To this end we follow |42] 
and introduce the weight 

w'd.//(^o,t) = e'^-(^("")'*)-^«) , (38) 

where x is the configuration space part of the phase space vector x. The diffusion constant 
may now be expressed in terms of the associated zeta function 1 42, ^] 

D = hm y{{x{xo),t) - xof) = hm 1(-^ + ^)(e'^-(^'(-o),t)-^o)) | 

The subscript diff will be suppressed form now on. 

We now consider the Lorentz gas obtained by infinitely unfolding the Sinai billiard. The 
relevant zeta function has a BER approximation (27) 

Zber{s) = 1 - (e-«^e"ift+-2/^2) ^ ^ _ (e"«^-l e^^^°^(^)d0),o, , (40) 

27r Jo 



where we assumed isotropy (which is not really necessary cf. refs. [p3| , p^), and A = yxf + x^. 
We only need to expand to second order in P 

Zber{s) = 1 - (e-^^) - ^(A^e-^) . . . (41) 
We can now evaluate the averages by means of the distribution pi3=o discussed in section Illb. 

^BEK(s) = l-(^ + 0(s2logs))-(^^^(logs + 0(l))j... . (42) 

keeping only the leading terms in R. Inserting tis expression into (^) we get as a result a 
diverging diffusion constant 

o = ,,. 1 r^- .''zMipm,, = + 0,1)) , (43) 
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in agreement with the suggested exact result |17|] and with numerical simulations |43]. This is 
one of the reasons why we expect that Z{s) = Zber{s) + 0{s^). 

This anomalous diffusion behaviour of this Lorentz gas is a consequence of the strongly 
intermittent properties of the Sinai billiard. Viewed as a Lorentz gas the diverging diffusion 
constant is due to the unbounded horizon, the point particle can make arbitrary long jumps 
between disks. 

d. A short cut to the trace formula 
In investigations of the accuracy of the trace formula for billiards the natural starting point 



is the boundary integral method |44, 29 1. Below we will, from this method, briefly derive the 
trace formula for the Sinai billiard. The reason for this exercise is that it will help us understand 
the interplay between neutral and unstable orbits in the next section. 

According to the boundary integral method the eigenvalues of the problem are those for 
which the following integral equation has a solution 

u{v{s)) = 2 / {T{s),T{s'))u{v{s'))ds' . (44) 

JS OUs 

The function u{r{s)) is related to the wave function according to 

«(r(.)) = ^^ ■ ^4^) 
The Green function is arbitrary but we now adopt the idea of ref. M] and use the one-disk 



Green function, the integral (44) should thus be performed only along the square boundary. 



The one disk Green function reads 

oo 

Giri,r2,Ae)='- J2 (hJ {kr^) + S,{kR)Hi {kr^)) Hj {kr2)e''^^'~^ , (46) 



where H^(z) are Hankel functions and ri, r2 and AO are explained in fig 3a. The phase shift 
function Si{kR) is defined by 

Using Poisson resummation we get 

oo 

G(ri,r2,A^)= ^ G^'^\n,r2, AO) , (48) 



m=— oo 



where 

G^'^Hn, r2, AO) ='-p {H^ikn) + S,{kR)H+{kri)) Hj {kr2)e''^^'+^^"'U^ . (49) 
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The standard semiclassical result are obtained by replacing all Hankel functions by their Debye 
approximation and taking the integral approximation of the sum (the higher terms in 

the Poisson resummation corresponds to creeping trajectories). The integral is then performed 
by stationary phase. 

We will now introduce an unnecessary approximation. The reason is just to make the 
following reasoning more transparent. We will assume that ri,r2 ^ R. This works well for 
small radii R because the end points in G{ri,r2, AO) is kept on the square boundary. In the 
semiclassical limit A: — > oo the Green function is 

j Gd{D) + Gd{Ri) V'iT^kR cos (pe-'^/'^Gd {R2 ) < arccos(i?/r) + arccos(i?/r') 
1 A0 > arccos(i?/r) + arccos(i?/r') ' 

expressing that in the lit region the Green function is a sum between the direct trajectory and 
one reflected on the disk. The lengths D, Ri and R2 are explained in fig 3a. Gd denote the 
asymptotic limit of the free (outgoing) Green function 

In ( |50[ ) we have neglected the transitional behaviour around AO ~ arccos(i?/ri) + arccos(i?/r2) 
which is referred to as the penumbra in this is the subject of section Ille. 
We can write eq (^) symbolically as the matrix equation 

(I-A)U = , (52) 

having a solution only when det(I — A) = 0. We can now express the integrated density of 
states as 

00 -I 
N{k) = ye{k- ki) = Im log det(I - A) (53) 

1 1 °° 1 

= Im tr log(I - A) = - -Im V -tr(A") , 

TT TT , n 

n=l 

where 

/f)G r)G r)G 

dsi... dsn^{r{s,), r{s2)) . . . -^(r(s„_i), r(.„))^(r(5„), r{s,)) . 

(54) 

Performing the normal derivative on (|50| ) we get (for large k) 



dG{ri r2, AO) ^ ^_ .^^ cos(gz3)Gd(I^) + (-ifc) cos{0R)Gd{Ri)V'^T^kRcos(j)e~'^/''Gd{R2) , (55) 
on 

provided of course that we are in the lit region AO < arccos(i?/ri) + arccos(i?/r2). The angles 
Od and Or are explained in fig 3b. 
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The integrals in (|54| ) will select classical orbits obeying the classical reflection law. After 



performing the last integral in (54) only the periodic orbits will remain. To compute the 
contribution from classical orbits we will need to perform the following integrals by stationary 
phase analysis around a stationary point s = sq 

J ds Gd{\v' - v{s)\)2{-ik) cos(e)Grf(|r(s) - r"|) ~ -Gd{\v' - r(so)| + |r(so) - r"|) • (56) 

Periodic orbits come in two different types. Orbits never bouncing on a disk are neutrally 
stable [neutral orbits for short) and enter in one-parameter families. Orbits bouncing at least 
once on a disk are unstable. In the small R limit the stability eigenvalues are approximately 

m 

rip r)J 

where li is the traveling length between the hits on the disk and (pi are the scattering angles 
for these hits. 

The contribution to tr(A")/n from a primitive unstable periodic p with length Lp hitting 
the square boundary Up times and traversed r = n/up times is 




4Mcos(0.),/^) e-'-. = . (58) 

The factor Up in front accounts for the number of periodic points on the square boundary and 
Xp counts the number of bounces along the orbit. The result above was obtained in the small 
R limit. For finite R the factor |Ap|^/^ should be replaced by |Ap|''/-^(l — l/Ap. The unstable 
orbits give rise to Gutzwiller's trace formula which can be recasted to the Gutzwiller-Voros 
zeta function. The neglect of the factor (1 — l/Ap is equivalent to the neglect of the higher m 
factors in (^) which is not our primary concern. 

A primitive neutral orbit (or rather a one-parameter family of neutral orbits) of length Lq 
and geometric width Dg and repetition number r = n/rig give the following contribution to 
tr(A'^)/n 

^ /2 1 \ „i(fcL,r-7r/4)+x,r7r 




4 V vr yfcV J ■ 
Now only n — 1 integrals are computed by stationary phase, the last integral gives rise to the 
factor Dg. 

The important lesson from this section is that the dichotomy between unstable and neutral 
orbits in the trace formula is a direct consequence of eq. (|50|). In the next section we will 
see how the neglect of the quantum fuzziness in the twilight zone between the reflection and 
direct region in the Green function ( |50| ) allows singularities to enter the Gutzwiller-Voros zeta 
function. 

e. How singularities sneaks into the Gutzwiller-Voros zeta function 
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We will consider scatterings in extreme forward angles so only the second term in the 
integral (|4|) (m=0) contributes. This means that we only need to evaluate the integral 

j POO 

G(°)(ri,r2,A0) = M Se{kR)H+{kn)H+{kr2)e'''^'de , (60) 

o J-OO 

where <C A9 < vr. We can still safely use the Debye approximation for the Hankel functions 
H^{kr) and H^ikr') 



y vrvz^ — P 

because the end points stay on the square boundary. However, the phase shift function needs 
a more careful analysis. The phase shift function Se{kR) is of unit modulus and we call the 
phase 'y{kR, £): 

5,(A;i?) = -:^^M^e^7(fci?/) . (62) 
^ ^ H+ikR) ^ ' 

The Green function now reads 

/•oo 

GW(ri,r2,A^) = / A(ri,r2)e^(^^^-*W)d^ , (63) 



where A(ri,r2), coming from the Hankel functions H^{kri) and {kr2), is a slowly varying 
amplitude so that the integral is in principle determined by the phase function 



qi{£) = -^[kR,l)-[^{krf -P-i arccos(£/(A;r))] - {kr'f -P-£ arccos(^/(fcr'))] +vr/2 , 

(64) 

and the stationary phase condition will simply read ^^{i) = A9 (subscripts denote differenti- 
ation). 

We will now investigate the phase of Si{kR) in detail. To this end we will use the uniform 
approximation for Hankel functions relating the phase of Hankel functions to the phase of Airy 
functions according to |45] 

S,[kR) = -Mz4±^5i(z4 ^ e2^(-)- = e^^(^'^^) (65) 
^ ' A\{-x) - iB\{-x) 

where 

J [^{^{kRY -e-£ arccos(£/(fci?)))]2/3 £ < kR 

'^'""''^ = j ■ log(^±4#^) - VF^W]^/^ l>kR ■ ^''^ 
If \kR/i — 1| <C 1 we can use 

^ = (2)i/3(^^_^) . (67) 

Using this expression for x one gets the so called transition region approximation. This ap- 
proximation fails to yield the Debye approximation as its asymptotic limit but the nice thing 
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is that there is a considerable overlap because whenever x <C ^^''^ we can use the transition 
region approximation and when x ^ 1 we can use Debye. (The relevant ^ for our consideration 
scales as £ ~ /c). 

6{x) is a complicated function but one has the following asymptotic expressions [45| 



3^/3 r(2/3)^ 35/3r(2/3)3 
^ ~ 87rr(l/3) 



e{x) = 7r/3 - ' Y"^^^ X - " ' + 0{x^) x^ibO (68) 



0(x) ~7r/2-^e-t(-^)'^'(l + 0(l/(-x)3/2)) x ^ -oo , (69) 

0(X) = ^ - + 0(x"3/2) X ^ +00 ( 70) 

so we obtain the Debye approximation when x — > oo as expected. 

In fig 4 we plot the function for some arbitrary choose values of ri/R = r2/R = 3 

and kR = 10 together with its Debye approximation. We have marked some critical values of 
^l{tj in the plot: 

• A^c; = arccos(i?/ri) + arccos(i?/r2) denotes the limit of the classically illuminated region. 

• A^rf is the limit of the standard semiclassical result on the direct side. It scales with kR 
as Md = AOci - 0{{kR)-^/^). 

• A6r is the corresponding limit on the reflection side. It scales as A6r = A6ci—0{{kR)^^/^). 

• The maximum of the curve occurs at A9„,ax = ^0d - 0{{kRy^/^ \og{kRf/^) and at £ = 
^■max- The second derivative of at the maximum is ^mi^max) = 0((A;i?)~'^/^(log kR)^/'^). 

If AO = AOmax + ^ and 6 is small, the Green function is still non zero but cannot be 
evaluated by stationary phase. This means that we cannot in an unambiguous way divide the 
Green function into one part associated with the direct ray and one with the reflected, as in 
(^0|), and something in the semiclassical description is definitely lost. 

Close to the maximum the result can be approximated by the integral 



oo 

e 

-oo 



i(5^+||*«.|(^-W)3)^^ ^ 27r|*OT/2|"i/3Ai(|^«^/2|-i/35) (71) 



Prom this we can get an estimate of the critical value of 5 when stationary phase approximation 
ceases to be valid. This happens when the argument of the Airy function is 0(1). We can 
define a critical 5 by 

5crit\'^m{^max)\=0{l) , (72) 

and we find that 5crit are of the order 

5,rit = 0{{kRr*l^{\ogkRfl^) . (73) 

We only expect the Airy function approximation to work well close to the maximum since the 
curve is very skew and badly approximated by a parabola. 
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In ref. |29] the authors choose to approximate the direct side i > k hy a Fresnel integral 
(this assuming constant slope of "^i) but the treatment on the reflections side is incomplete. 

Let us now consider the family g„ of periodic orbits indicated in fig 1. We see that if n is 
big enough the unstable periodic orbit will interfere with the neutral orbit and for large enough 
n the corresponding terms in the Gutzwiller trace formula are suppressed. More precisely, it is 
the stationary points closest to the disk (cf. the one marked with a cross in fig 1) that suff'ers 
most of interference with the edge of the corresponding family of neutral orbits. In the same 
way as there is no longer an unambiguous division between a direct and a reflected ray there 
is no longer a corresponding division between neutral and unstable orbits in the trace formula, 
cf. section Illd. It is the neglect of this effect that allows the branch point singularity in the 
Gutz wilier- Voros zeta function. 

The critical n = ricrit when this interference problem gets crucial scales as 

n,Ht ~ (fei?)'/' . (74) 



which is the same as the critical threshold of ref. |29| above which the majority of periodic 
orbits suffer from 'diffraction' effects. 

It is this type of entanglement which makes semiclassical quantization of systems with 
mixed phase space so difficult. Our findings are very similar to those of ref.[^] concerning the 
Stadium billiard. 



I would like to thank Karl-Erik Thylwe for discussions on circle Green functions. To result 
in section Ille will appear in a joint publication. This work was supported by the Swedish 
Natural Science Research Council (NFR) under contract no. F-FU 06420-303. 
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Figure captions 



Figure 1: The Sinai billiard (upper left corner) is the unit cell of the associated Lorentz gas. 
Two transparent direction are indicated (shaded areas). The family g„ of periodic orbits accu- 
mulating towards the horizontal transparent direction is also indicated. 

Figure 2: The intermittent map discussed in section Illb. The sequence {cn} is are defined as 
the inverse images of the critical point as indicated. 

Figure 3: The semiclassical limit of the Green function between two points on the border can 
be divied into two parts. One associated with the direct trajectory and one that is reflected on 
the disk. 

Figure 4: The function ^'^ discussed in section Ille together with its Debye approximation. 
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